Capture numbers and island size distributions 
in models of submonolayer surface growth 
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The capture numbers entering the rate equations (RE) for submonolayer film growth are de- 
termined from extensive kinetic Monte Carlo (KMC) simulations for simple representative growth 
models yielding point, compact, and fractal island morphologies. The full dependence of the capture 
numbers a s (0,T) on island size s, and on both the coverage O and the V = D/F ratio between 
the adatom diffusion coefficient D and deposition rate F is determined. Based on this information, 
the RE are solved to give the RE island size distribution (RE-ISD), as quantified by the number 
n s (O, r) of islands of size s per unit area. The RE-ISDs are shown to agree well with the correspond- 
ing KMC-ISDs for all island morphologies. For compact morphologies, however, this agreement is 
only present for coverages smaller than ~ 5% due to a significantly increased coalescence rate 
compared to fractal morphologies. As found earlier, the scaled KMC-ISDs n s s 2 /0 as a function of 
scaled island size x — s/s approach, for fixed O, a limiting curve foo(x, O) for T — > oo. Our findings 
provide evidence that the limiting curve is independent of O for point islands, while the results for 
compact and fractal island morphologies indicate a dependence on 6. 

PACS numbers: 81.15.Aa,68.55.A-,68.55.-a 



I. INTRODUCTION 

The kinetics of submonolayer nucleation and island 
growth during the initial stage of epitaxial thin film 
growth has been studied intensively both experimentally 
and theoretically for more than three decades (for re- 
views, see Refs. U0, and references therein). Impor- 
tant aspects of the growth kinetics in the submonolayer 
growth regime can be described by the rate equations 
(RE) approach^ This approach has proven to be very 
valuable in inorganic thin film growth. Interestingly, 
many of the theoretical concepts developed for thin film 
growth kinetics of inorganic materials, have shown re- 
cently to be very valuable also for applications in organic 
thin film growth^— This is due to the fact that these 
concepts often are not specifically referring to particular 
materials. Instead, they take into account the key mech- 
anisms involved in the complex interplay of deposition, 
evaporation, diffusion, aggregation and dissociation from 
a general viewpoint. 

Parameters entering the RE are the capture numbers 
tr s (9, r), which describe the strength of islands of size s 
to capture adatoms at a coverage 9 and ratio D/F = T 
of the adatom diffusion coefficient D and deposition flux 
F. The dependence of the capture numbers on s has been 
studied for various T but only for one or a few 9 values. 
In this work we present a systematic study of the full 
dependence on both 9 and V for different types of island 
morphologies and the case, where detachments of atoms 
from islands can be neglected, corresponding to a critical 
nucleus of size i = 1 . This is motivated by the following 
questions, which have not been thoroughly answered yet: 

(1) If the cr s (9,r) are known, do the RE then predict 
correctly the number density n s (9,r) of islands 
of size s, that means the island size distribution 
(ISD)? This question indeed was earlier posed by 



Ratsch and Venables 1 -^ as well as Evans et al& The 
answer to this question is not obvious, since the 
RE with known capture numbers a s (Q 1 T) neglect 
many-particle correlation effects^ spatial fluctua- 
tions in shapes and capture zones of islands, and 
coalescence events that, despite rare in the early- 
stage growth, can have a significant influenced 

(2) Is there a simple functional form of the a s (9, T), in 
particular, is there a scaling of these capture num- 
bers with respect to an effective capture length as 
suggested by a self-consistent treatmen t 14 ' 15 based 
on the RE? Do the <r s , when scaled with respect 
to their mean <r, depend for large T on the scaled 
island size s/s only, as suggested by Bartelt and 
Evans±2? 

In previous studies it has been found that the scaled 
ISD s 2 n s /Q as a function of scaled island size x = s/s 
approaches, for fixed coverage 9, a scaling function f{x) 
for large T. Early simulations suggested that f(x) is in- 
dependent of 9 and moreover not sensitive to the island 
morphology. However, later results showed^ 1 ^ that the 
morphology has an influence on the form of f(x). In fact, 
one would expect the scaling function f(x) to become in- 
dependent of 9 if the RE with known capture numbers 
correctly predict the ISD, and if the scaled capture num- 
bers a /a as a function of s/s become independent of 9 
for large T. Under these assumptions, an explicit relation 
was proposed by Bartelt and Evans) 16 ' 17 which connects 
the scaling function of the capture numbers with the scal- 
ing function of the ISD. We hence address the following 
further questions: 

(3) Is the scaling function f(x) independent of 9 for 
large T? What is the influence of the island mor- 
phology? Can the relation between the scaled ISD 
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f(x) and the scaling function for the capture num- 
bers be confirmed? 

The RE treatment is based on a coupled set of sim- 
ple rate equations describing the time evolution of the 
adatom density tlx and the number density n s of islands 
with size s > 2, if spatial correlations among islands dur- 
ing growth are neglected. Taking into account direct im- 
pingement of arriving atoms at the border of islands, the 
RE for the case i — 1 read 

4~jT" =(! - ) - ZFoin\ - Tni ^ a s n s 

S>1 



2kiUi 



S>1 



1 dn s 
F~df 



-Tni (cr s _in s _i - cr s n s ) 



K s _in s _i 



2,3,.. 
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(2) 



These equations refer to the pre-coalescence regime where 
only adatoms are mobile and it is presumed that re- 
evaporation of atoms and atom movements between the 
first and second layer can be disregarded. Moreover, 
adatoms arriving on top of an island are not counted, 
i.e. s in a strict sense refers to the number of substrate 
sites covered by an island (or the island area). The cov- 
erage 8 entering Eq. {jXJ) is given by = J2 s >i sris = 
1 — exp(— Ft) and takes into account that adatoms are 
generated by deposition into the uncovered substrate 
area (as common in the literature in this field, we set 
the length unit equal to the the size of the substrate lat- 
tice unit). The terms with 01 (0) and ki(0) describe the 
nucleation of dimers due to attachment of two adatoms 
by diffusion and due to direct impingement, respectively. 
The term oc nia s n s describes the attachment of adatoms 
to islands of size s > 1, and the term oc K s n s the direct 
impingement of deposited atoms to boundaries of islands 
with size s. For the idealized point island model, s refers 
to the total number of atoms that arrived at a point, and 
(1 — 0) in Eq. ([T]) is replaced by one (no covered sub- 
strate area) . For a unified discussion of capture numbers 
and the ISD we formally set = Ft for the point island 
model. 

Introducing the total number density of stable islands 
N and the average capture number a, 



S>1 



(3) 



S>1 



a reduced set of equations for Ui(0) and iV(0) can be 
derived from Eqs. (HJ and ^ within the RE treatment. 
These equations predict the scaling relation N oc T~ x 
with the scaling exponent \ — 1/3 . 18 i 19 This relation 
has been successfully validated by several growth exper- 
iments in the past and applied to extract adatom diffu- 
sion barriers and binding energies in metal epitaxy. A 
discussion of many of these experiments can be found in 
Ref. l2|. Recently, the relation has also been applied in or- 
ganic thin film growth &± An extended RE approach for 



multicomponent adsorbates^S was recently suggested to 
determine binding energies between unlike atoms from is- 
land density data.— 

More detailed information on the growth kinetics is 
contained in the ISD. If the full dependence of the ISD 
n s (0, r) on and T is mediated by the mean island size 
s(0,T), the ISD should obey the following scaling form, 
as first suggested by Vicsek and Family^ 



2 (e,r) 




n s (6,r) = / 



»(e,r) 



(4) 



Here the scaling function f(x) must fulfill the conditions 
J Q f{x)dx = J Q xf(x)dx — 1. The scaling behavior 
was found to give a good effective description for large 
r. More precisely, the curves s 2 n s /0 as a function of 
x — s/s approach a limiting curved 



lim —n xS (Q,T) = f 00 (x,Q). 
r— >oo t? 



(5) 



Previous studies for a few fixed values suggest that 
foo{x,Q) is independent of 0. 

An explicit expression for the scaling function f(x) 
with shape independent of was suggested by Amar 
and Family^ 



f(x) = CiX 1 exp (—idiX 1 ^^ . 



(6) 



The parameters entering this scaling function depend on 
the size of the critical nucleus i, which allows one to deter- 



mine i in experiments 
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Equation (6J was believed 



to be independent even of the morphology 24 , but this has 
later been questionedj^^ 

Based on a continuum limit of the RE @ and scaling 
assumptions for the capture numbers and a neglect of 
the ©-dependence, an expression for the limiting curve 
foo{x) was derived by Bartelt and Evans ; 16 ' 17 



foo(x) = /oo(0) exp \ J dy 



(2 Zoo -l)-C t ' ot (y) 
Ctot(y) - zy 



, (7) 



where = 9(lns)/9(ln0), and C to t(x) is a linear com- 
bination of the scaled capture numbers Coo(x) — a s /a 
and scaled direct capture areas K^^x) = k s /k. The co- 
subscript indicates that the large T — > 00 limit should 
be taken. As pointed out by Bartelt and Evans, Ct t{x) 
should be well approximated by the scaled capture num- 
bers alone, C to t{x) w (1 — 0)6*00 (x). In Appendix [Al we 
show that in fact it holds Ct t(x) ps Coo(x). The two 
conditions for foo(x) (normalization and first moment 
equal to one) imply that Coo(0) = (1 — z)//oo(0) and 

j™dxc 00 (x)f 0? (x) = im 

It is interesting to note that a semi-empirical form, 
which has a structure similar to Eq. © has been sug- 
gested recently by Pimpinclli and Einstein^ 7 - for the dis- 
tribution of capture zones A as identified by Voronoi tes- 
sellation, 



Pf) = Cf3 a exp(— dpa 2 ) 



(8) 
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where a = A/ A is the rescaled capture zone with respect 
to the mean A and (3 = i + 1 (see also Ref. (28|). This 
distribution corresponds to a generalized Wigner surmise 
from random matrix theory. The parameter ft = i + 
1 and the functional form, however, are controversially 
discussed. 29 - 30 



Besides this recent progress in predicting functional 
forms of capture zone distributions, there are only a few 
studies so fa r 13 ' 31 ' 32 that address the problem whether 
an integration of the RE (PJ and ([2]) can yield correctly 
the ISD for different cluster morphologies in the pre- 
coalescence regime. For an integration of the RE a re- 
liable determination of a s (@,T) is needed. Four gen- 
eral approaches have been followed for this purpose: (i) 
Within a self-consistent ansatz one can solve the diffusion 
field around an island and derive determining equations 
for the capture numbers by equating the attachment cur- 
rents of the diffusion field and the REi 14 i 15 (ii) By mod- 
eling the island growth with the level set method^ 3 - one 
can analogously equate attachments currents and deter- 
mine the capture numbersi 31 ' 32 (iii) Balancing the de- 
position rate FA s n s into the mean capture zone A s of 
islands of size s with the RE expression Da s ri\n s for the 
attachment rate to these islands, yields a s ~ A s /Tn\. 
This means that the capture numbers can be approxi- 
mately calculated from a determination of the A s , e.g. by 
Voronoi tessellation! 17 i 34 ~ aj£ (iv) In simulations, where the 
individual attachments are followed, the capture num- 
bers can be calculated from the mean number of attach- 
ments M s to island of size s during a time interval At 
[sec Eq. fljJJ) and the discussion in Sec. Hi]. 16 



The paper is organized as follows. First we describe 
in Sec. [TT] the method used to generate point, compact 
and fractal island morphologies, and the method for de- 
termining the capture numbers as function of island size 
and coverage. In Sec. Mil we discuss the results for the 
capture numbers and compare these with the prediction 
of the self-consistent theory. In Sec. IIVI we analyze the 
mean island and adatom densities for the different is- 
land morphologies and discuss their prediction by the 
self-consistent RE and the RE based on the capture num- 
bers determined in the KMC simulations. In Sec. [V] we 
demonstrate that the ISD is successfully predicted by the 
RE as long as coalescence events can be neglected. These 
coalescence events are relevant already for small cover- 
ages > 0.05 for compact morphologies, while they turn 
out to be much less important for fractal morphologies. 
The reason for these differences are reduced coalescence 
rates for fractal island morphologies because of a screen- 
ing effect jii^i Finally, we study in Sec. lVII the behavior of 
the scaled capture numbers and scaled ISDs in the limit 

r -> oo. 



II. SUBMONOLAYER GROWTH: MODELS, 
MORPHOLOGIES AND SIMULATIONS 

The KMC simulations are performed with a first reac- 
tion time Monte Carlo algorithm 38,39 on a square lattice 
with L x L sites. In this algorithm, two times tf and td 
are randomly generated from the exponential probability 
density ip(r) = 7exp(— "yr), where 7 = L 2 F for rp, cor- 
responding to a deposition process, and 7 = ADL 2 n\ for 
td, corresponding to one of the possible diffusive jumps 
of adatoms. If tjj < tf, the simulation time is incre- 
mented by td and one of the L 2 n\ adatoms is selected 
randomly and moved to a randomly selected vacant near- 
est neighbor site. If Tp < td, the simulation time is 
incremented by Tp and one of the L 2 sites is randomly 
chosen. If this site is vacant, an additional adatom is 
deposited on this site, while, if the site is occupied, no 
deposition takes place. 

With respect to the formation of islands we consider 
three simple growth models that are representative for 
the different types of island morphologies in the case of 
i = 1. Fractal islands are generated by applying "hit 
and stick" aggregation, that means an adatom having 
another atom as nearest neighbor becomes immobilized. 
Compact island morphologies are produced by lettin g is- 
lands grow spirally into a quadratic form as in Ref. |40j, 
meaning that each adatom attaching to an island is dis- 
placed to the corresponding tip of the spiral. Point island 
morphologies are generated by displacing an adatom at- 
taching to an island to the site representing the island, 
while bookkeeping the total number of aggregated atoms 
for the island size. 

To calculate the capture numbers a s at the cover- 
age 0, we use the following procedure which is based 
on the method outlined in Ref. Il6l : Each simulation 
run is stopped at coverage and the number densities 
n s = N s /L 2 , s — 1, 2, . . . are determined, where N s are 
the numbers of monomers (s = 1) and islands (s > 1). 
Then the simulation is continued for a time interval At 
without deposition and the following additional rules are 
implemented: (i) if an adatom is attaching to an island of 
size s > 1, a counter M s is incremented and the adatom 
thereafter repositioned at a randomly selected site on 
the free substrate area (i.e. a site which is neither cov- 
ered nor a nearest neighbor of a covered site); (ii) if two 
adatoms form a dimer, a counter Mi is incremented and 
the two adatoms thereafter are repositioned randomly as 
described in (i). In this way a stationary state is main- 
tained at the coverage 0. The mean attachment rate per 
unit area to islands of size s is M s /(L 2 At), and equat- 
ing this with the expression Da s n a n\ from the RE (|1I2[) 
yields 



L 2 AtD 



n s ni 



1,2,. 



(9) 



Averaging over many simulation runs (configurations) 
gives <T S (0, r). The K s £1X6 determined from the lengths 
of the islands boundaries, which are simultaneously mon- 
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FIG. 1. Capture numbers as a function of s for F = 10 7 and four different fixed coverages for the models representing (a) point, 
(b) compact, and (c) fractal island morphologies. The filled symbols refer to the a s obtained from the KMC simulations and 
the open symbols to the results of the self-consistent theory according to Eqs. (|10[1 (with the R 3 taken from the simulations, 
see text and Fig. [3|. 



(a) point islands (b) compact islands (c) fractal islands 








FIG. 2. Mean capture number a and mean island size s as a function of 6 for the four simulated V values and the models 
representing point, compact and fractal island morphologies. 



itored during the simulation and averaged for each size 
s. 

The continuous-time Monte Carlo (KMC) simulations 
are performed on a square lattice with periodic bound- 
ary conditions and L x L = 8000 x 8000 sites for four 
different T = 10 5 , 10 6 , 10 7 , and 10 s . For each value of 
an average over 10 s nucleation/attachment events was 
performed. 



III. CAPTURE NUMBERS 

The direct capture areas for point islands on a square 
lattice are given by n s = 4. For compact and fractal 
islands the k s increase as ~ y/s and ~ s, respectively, 
and their dependence on and T is very weak. 

Representative results for the capture numbers a s are 
shown in Fig. Q] as a function of s for fixed T — 10 7 and 
four different coverages for the (a) point, (b) compact, 
and (c) fractal island morphologies. For the other sim- 
ulated r values, a similar behavior was obtained. The 
mean tr(0, V) [see Eq. ©] as a function of for all sim- 
ulated r values is displayed in Fig. [21 together with the 



mean island size s(0,r). These functions are later used 
in Sec. IVII when investigating the scaled capture numbers 
a s in connection with the scaled island densities in the 
limit r — > oo. 

A common feature for all morphologies in Fig. [T] is a 
linear increase of a s with s for large s > s. It can be 
understood^ from the proportionality of the a s to the 
mean capture zone areas A s , and the fact that large is- 
lands typically exhibit large A s , which led to the stronger 
growth of these islands. Since a twice as large capture 
zone gives on average rise to a twice as large island, it 
holds A s ~ s and hence a s oc A s ~ 

With respect to the dependence on the coverage 0, the 
cr s in Fig. [T]have a quite different behavior for the three 
morphologies in the regime s > s: While for the point 
islands the a s decrease with 0, they are almost inde- 
pendent of for the compact islands, and they increase 
with for the fractal islands. Main reason for these dif- 
ferences is that for point islands the number density N 
continues to increase with (that means time t = Q/F) 
due to ongoing nucleation of new islands, while for com- 
pact and fractal morphologies, N tends to saturate for 
larger 0, with less pronounced saturation in the com- 
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FIG. 3. Mean radii of gyration R s of islands of size s for the 
models representing compact and fractal island morphologies. 
The straight lines in the double-logarithmic plot indicate the 
power-law behavior for large s. 



pact case (see Sec. II VI below). During the growth in the 
point island model, a large capture zone surrounding a 
large island is, compared to the other two morphologies, 
more frequently destroyed by a nucleation event in this 
zone, and the A s oc er s thus decrease with for fixed 
s > s. Due to the higher nucleation rate and the missing 
spatial extension of islands in the point island model, the 
corresponding a s are much smaller than for the compact 
and fractal morphologies. The larger island extension 
and the strong capture of adatoms by finger tips in the 
case of fractal islands lead to about five times larger <j s 
in comparison to the compact islands. 

The differences with respect to the dependence are 
also reflected in the behavior of cr(0, T) in Fig. [21 In fact, 
when considering the scaled capture numbers a s /a, the 
dependence for s > s becomes qualitatively the same 
for all morphologies (increase of a s /a with 0, see Sec.lVll 
below). For small s < s, the curves in Fig. [T] show a non- 
linear dependence of <j s on s for all morphologies i 35 ' 41 ! 42 
By combining the linear function for large s with a poly- 
nomial at small s, we fitted the results for a s for all sim- 
ulated and r values. These fits, together with corre- 
sponding fits for the k s , were used to integrate the RE 
© and ©. 

The mean island size s in Fig. [5] reproduces the be- 
havior seen in many earlier studies.— In the point island 
model the straight lines in the double logarithmic repre- 
sentation are in agreement with s ~ Z with z = 2/3 as 
predicted by a scaling analysis of the reduced RE ^ 43 i 44 
In the case of the compact and fractal island morpholo- 
gies, the slope z(Q) — <91ns(0,r)/<91n0 increases with 
and approaches z ~ 1 for both island morphologies. 
This is consistent with a saturation (©-independence) 
of the island density for large in the pre-coalescence 
regime, s ~ Q/N ~ 0. 

In the self-consistent theory-^ 4 - the capture numbers are 




FIG. 4. Scaling plot of (f — <d)a s as a function of R a /£ for the 
model representing compact island morphologies, with both 
R B and £~ 2 = 2aini + Yls>2 o~ s n B determined from the KMC 
simulations, (a) for various V and fixed O = 0.2, and (b) for 
various Q and fixed F = 10 8 . The insets in (a) and (b) show 
the corresponding results for the model representing fractal 
island morphologies. The solid lines represent the specific 
functional form in Eq. (|f0l) predicted by the self-consistent 
theory. 



given by 



sc = 2tt R s Ki (R s /£) 

° s (i - 0) e K (i? s /o ' 
r 2 = 2^ + 5>> s , 

s>2 



(10a) 
(10b) 



where R s is the effective radius of an island of size s, 
Kq and K\ are the modified Bessel functions of order 
zero and one, respectively, and £ is the adatom capture 
length (mean linear size of depletion zone around an is- 
land). The factor (1 — 0) in Eq. (fTUl) . which was not 
given in the original derivation in Ref. 14j arises from the 
fact that the adatom current to a (circular) island of size 
s is 2-kR s D d r n\{r)\ r= R g , where h\(r) = m(r)/(l — 0) 
is the adatom density with respect to the free (uncov- 
ered) surface area, and ni(r) is the local form corre- 
sponding to the global mean value m appearing in the 
RE [see also Ref. |45| for the additional factor (1 — 0)]. 
For R s < £, a s ~ 2n/[{l - 0) ln(£/i? s )], and for R s » £, 

a s ~ 2tt/(1 - Q)(R./0- 

To determine the al c , the RE |T]) and (J2J) are numeri- 
cally solved with initial conditions n s — at time t = 
and a cutoff value s c so that n s can be safely neglected 
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FIG. 5. KMC results (symbols) for the adatom density n\ and island density N as function of G at F — 10 r for (a) point, 
(b) compact, and (c) fractal island morphologies in comparison with the RE solutions, when using the direct capture areas 
K a and capture numbers a 3 from Sec. Mil (solid lines) and when using the capture numbers al c from the self-consistent theory 
according to Eq. (flO|) (dashed lines) . 



for s > s c . In each integration step the implicit Eq. (fTU)) 
is solved for the a^ c . The results become sensitive to the 
island morphology via the dependence of R s on s in this 
approach. For point islands we take R s = 1 correspond- 
ing to one lattice constant. For the compact and fractal 
island morphologies, we determined the mean radius of 
gyration of islands of size s, as shown in Fig. |3l The 
straight lines in the double-logarithmic representation 
give R s ~ 0.42s 1/2 (compact islands) and R s ~ 0.47s 57 
(fractal morphologies) for large s. To compare the <r^ c 
with the a s obtained from the KMC simulations, we used 
the full dependence of the R s on s, i.e. including the small 
s behavior, in our integration of the RE. The results from 
the self-consistent theory are shown in Fig. [T] (open sym- 
bols). As can be seen from the figure, the af 1 deviate 
strongly from the KMC results, both in their size and 
in their functional form. In particular the self-consistent 
theory underestimates the capture numbers for large s, 
as known from earlier work in the literature^. 

It is interesting to see, whether the scaling of (1 — 
0)cr s (0,r) with R s /£ is valid, if the a s and n s from the 
KMC simulations are used in the expression for £ -2 in 
Eq. ([TOb]) . In this case the linearization step used in 
this theory for deriving a linear diffusion equation for 
the local adatom density n\{r) could be reasoned, i.e. 
the step, where the term 2uini(r) 2 +X) s >i a sni(r)n s (r) 
is replaced by ni(r)£ -2 with £~ 2 = 2aini + J2s>i a « n *_ 
given by the mean (r-independent) densities (see Ref. [lj 
for details). In Fig. - @)c s (0,r) is plotted as func- 
tion of R s /£ for the models representing compact and 
fractal island morphologies. Figure IDJa) shows that in- 
deed a data collapse is obtained for different T values 
at fixed 6. However, with respect to the 0-dependence, 
tested in Fig. HJb), no scaling behavior is found. This 
indicates that the linearization step in the self-consistent 
theory leads to the unsatisfactory capture numbers. It 
has been shown that correlation effects between island 



sizes and capture areas need to be taken into account 
to improve theories for capture numbers and island size 
distributions. This can been achieved by considering the 
joint probability of island size and capture area i 42 i 45 ~ — 



IV. ADATOM AND ISLAND DENSITIES 

Numerical integration of the RE with the k s and a s 
from Sec. Mil gives an excellent description of the adatom 
density ni and of the island density N as a function of O 
and r for all island morphologies in the pre-coalescence 
regime. This is demonstrated in Fig. [5j where ni and N 
from the KMC simulation (open squares) and RE solu- 
tion (solid lines) are plotted as function of 9 for T = 10 7 . 
For compact and fractal island morphologies the KMC 
data for N steeply fall for coverages larger than 15% 
(compact islands) and 30% (fractal islands) because of 
island coalescences. Small deviations of the RE solution 
for m can be seen close to its maximum, where it slightly 
underestimates the adatom density. The agreement for 
the other simulated T values is of the same quality. As 
known from previous studies^ i n ' 14 i 45 the RE predict n\ 
and N quite well also, when using the self-consistent cap- 
ture numbers from Eq. (I10[) . The corresponding solutions 
are drawn as dashed lines in Fig. [5] In view of the discrep- 
ancies discussed in Sec. IIIII this good predictive power of 
the RE under use of the self-consistent capture numbers 
ct^ c is surprising. 



V. ISLAND SIZE DISTRIBUTIONS 

Since the a s s c (Q, T) deviate strongly from the a s (Q, T), 
the RE with self-consistent capture numbers fail to pre- 
dict the ISD. This failure was reported already when the 
self-consisting theory was developed*^ In the following 
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FIG. 6. Simulated island size distribution (circles) in com- 
parison with the RE solution (lines) for three different G at 
r = I0 7 for the models representing (a) point and (b) fractal 
island morphologies. 



FIG. 7. (a) ISD and (b) sub-ISD for compact islands obtained 
from the KMC simulations (circles) in comparison with the 
RE solution (lines) for T — 10 7 and three different coverages 
6. 



we therefore do no longer consider the self-consistent the- 
ory, but concentrate on the principal questions whether 
the RE with the capture numbers cr(0,r) are success- 
ful in predicting the ISD, and if so, whether in the limit 
r — > co the asymptotic form ([7]) for the scaling function 
becomes valid. In this section we address the first of these 
two questions. 

Representative results for the ISD (symbols) in com- 
parison with the RE predictions are shown in Fig. [6] for 
r = 10 7 and three different coverages 0, for point and 
fractal island morphologies. The excellent agreement be- 
tween the RE predictions and the KMC data in that 
figure is also found for the other simulated V values. As 
was shown in Ref. for the fractal island morpholo- 
gies, a x 2 test with a standard significance level of 5% is 
passed up to a coverage of 6 = 0.18. For larger 0, coales- 
cence events, not included in the RE approach, become 
relevant. 

For the compact island morphologies, a good agree- 
ment of the KMC data with the RE prediction is obtained 
up to coverages of about = 0.05 only, see Fig. [7k. The 
reason for the discrepancies are coalescence events that 
become important already for small > 0.05, in con- 
trast to what one may conclude from the behavior of the 
mean island density shown in Fig. [SJ where coalescences 
seem to be irrelevant up to coverages of about 15%. One 



can take out the coalescence effect in the calculation of 
the ISD by following the islands in the simulations and 
by counting coalesced islands as if they were separated. 
The islands identified in this way were referred to as sub- 
islands and the resulting ISD as sub-ISD in Refs. |H and 
Eol . In the same way as described in Sec. El we deter- 
mined the k' s and a' s for the sub-islands and integrated 
the RE ([1) and © with these input quantities. As shown 
in Fig.[jjb), these RE results for the sub-ISD give again 
excellent agreement with the KMC data. 

That coalescence events are much more frequent for 
compact than for fractal islands is shown in Fig. [8j where 
we plotted the fraction of the coalesced islands as a func- 
tion of both for the compact and fractal island mor- 
phologies. This fraction was determined by dividing the 
total number of coalescences up to the coverage by the 
total number of islands at this value, i.e. islands which 
have undergone more than one coalescence are counted 
with their corresponding multiplicities. As can be seen 
from Fig. the fraction of coalesced islands for compact 
islands has already at = 0.05 reached a level compara- 
ble to that found for the fractal islands at = 0.2. 

The reason for the less frequent coalescences of fractal 
islands is that two approaching fractal islands can avoid 
each other for some time, because fingers of one islands 
grow into breaches between fingers of the other island. 



8 



10" 



:io' 



LO 



■ compact islands 
» fractal islands 




10 



FIG. 8. Fraction of coalesced islands as a function of the 
coverage <d for T — 10 7 . 



When a finger enters a breach, its further growth slows 
down because of the shielding inside the breach. This 
screening effect and its consequence for coalescences has 
been discussed earlier in the literature.—^ A quantita- 
tive analysis of the coalescence behavior of compact and 
fractal island morphologies is given in Appendix [B] 



VI. LIMITING BEHAVIOR FOR r — s- oo 

Based on our first key finding that for all morpholo- 
gies and for all coverages in the pre-coalescence regime, 
the ISDs from the KMC simulations are successfully pre- 
dicted by the RE we now turn to the question, whether 
the scaled ISDs approach the asymptotic form ([7]) in the 
limit r — y oo . 

To answer this question is not easy because of 
various subtleties, which let us revisit the deriva- 
tion of the scaling function by Bartelt and Evan o 16 i 17 
in Appendix [X] As mentioned in the Introduction, 
Eq. ([4]) for the limiting curve of the scaled ISD foo(x) 
should be valid if foo(x,Q) from Eq. ([5]) is inde- 
pendent of 0. This is the case if Coo(a;,8) = 
limr->oo 0xs{Q, r)/<r(0, T) for the scaled capture num- 
bers and z ao (Q) — limp-j-oo 91ns(0, r)/91n0 also have 
©-independent limits. A further requirement for the va- 
lidity of Eq. is that lim r ^ooK(0,r)/s(0,r) = 0, 
where R = N^ 1 X) s >i K s n s is the mean direct capture 
area. This condition can be expected to be fulfilled for 
compact and point island morphologies and is in fact 
the reason, why the scaling function of the direct cap- 
ture areas should not enter the RE prediction ([7]). If 
limr-s-oo k(0, T)/s(0, T) > 0, foo{x,Q) can be expected 
to depend on and one would need to solve the semi- 
linear partial differential equation (|A6|) for fco(x,Q). 
Note that the n s cannot increase stronger than linearly 
with s, and accordingly R should not increase more than 
linearly with s. 

In interpreting numerical results for finite T, we have 
to pay attention to the fact that for smaller larger T 
values are needed to approach the limiting curves. This 



10 




10' 



°r=io 0=0.2 



(a) point islands 

0=0.02 
0=0.05 
0=0.10 
0=0.20 

r=io s 



2 

_J L 

1.5 



(b) fractal islands 

° 0=0.02 
o 0=0.05 
* 0=0.10 
o 0=0.20 




FIG. 9. Scaled island size distributions n s s 2 /Q as function 
of the scaled island size x = s/s for fixed V = 10 s and four 
coverages Q, for (a) point and (b) fractal island morphologies. 
The insets show the scaled ISDs at fixed <d = 0.2 and the four 
simulated F values. The line in (a) is a fit to the data for 
9 = 0.2 and T = 10 s and agrees with the analytical result 
Eq. 1|7|). when using the line in Fig. HOf al as the estimate for 



is because s(0,T) must become large enough to reach 
the "continuum limit" (and larger T are needed to ob- 
tain the same s at smaller 0), and because the relation 
ni ~ (1 — Q)/TaN, used in the derivation of Eq. ((JJ, 
should be obeyed. This relation is usually referred to 
as the quasi-stationary condition, since it follows from 
balancing the adatom attachment rate DanN to islands 
with the deposition rate ^(1—0). However, as was shown 
earlier ^ the relation is also valid for small values in the 
regimes, where relative changes of N are still large and 
have not leveled off. A refined scaling analysis^ yields 
that, for i = 1 as relevant here, the relation holds for 
(02r)V3 > i ; implying again that for smaller larger 
r are needed to identify the limiting behavior. 

Figure [TJ] shows n s s 2 /<d as a function of s/s for T = 10 8 
at four different coverages for the (a) point and (b) frac- 
tal island morphologies. In the insets the scaled ISDs are 
shown for a fixed coverage = 0.2 and different T. In 
the case of the point island morphologies, the data sug- 
gest the existence of a ©-independent limiting curve, in 
agreement with previous findings^ For the fractal island 
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FIG. 10. Scaled capture numbers cr s /<7 as function of the 
scaled island size s/s for fixed F = 10 s and four coverages 0, 
for (a) point and (b) fractal island morphologies. The insets 
show the scaled capture numbers at fixed <d = 0.2 and the 
four simulated F values. The line in (a) marks the solution 
obtained from a numerical integration of Eq. pip when using 
the fit (line) to the scaled ISDs curve for F = 10 8 and 6 = 0.2 
in Fig. Ha). 



morphologies, the scaled ISD for different 9 show no clear 
signature of a 0-independent limiting curve. Based on 
the tendency of the simulated data for different and T 
to become slightly closer to each other for larger and 
F, one may conjecture that also in this case a limiting 
curve would be reached at larger T values. However, the 
fact that for each fixed 0, the curves at large T are al- 
most overlapping suggests that these are good estimates 
°f foo( x i®)' O ur conclusion is therefore that it is not 
likely that a ©-independent limiting curve exists for the 
hit-and-stick model used here for the fractal island mor- 
phologies. 

This conclusion is further corroborated by the fact 
that the scaled direct capture areas exhibit a nearly lin- 
ear dependence on x for the fractal islands (not shown). 
Thus we encounter the case here, where the scaled di- 
rect capture areas k(0, r)/s(0, T) appear to approach 
a non- vanishing limit for r — > oo, which would mean 
that in a strict treatment, Eq. (J7J can no longer be ap- 
plied. If one considers foo(x,Q) to depend only very 
weakly on ©, df OD (x,0)/dO ~ 0, we could replace 



Coo(z) by C tot (x) = (a) + p^QK^x, ©)/(l - 0), 
where K OQ (x 1 ©) = limp— yoo k>xs 

(0,r)/S(0,r) is the lim- 
iting curve for the scaled direct capture areas and poo — 

iim r ^ooK(0,r)/s(0,r) > o. 

Our conclusions drawn with respect to the scaled ISDs 
of the point islands are consistent with the behavior of the 
scaled capture numbers, which are shown in Fig. llOf a) for 
the same and T values as in Fig.^a). In this Fig.QIJa) 
an approach to a ©-independent limiting curve C oc (x) 
can be seen. In the case of the fractal island morphologies 
by contrast, an approach to a ©-independent limiting 
curve cannot be clearly identified, which gives further 
evidence that the foo{x, ©) are dependent on 0. 

In order to test the validity of Eq. ([7J for the point 
islands, we set = 2/3 (see Sec. IIII[) and used a fit to 
the scaled ISD for T = 10 8 and = 0.2 in Fig.^a) as an 
estimate for foo(x). The fit, which fulfills the constraints 
of normalization and normalized first moment, is shown 
as line in this figure. We then estimated Coo(x) based 
on this fit by rewriting Eq. (|A6|) from Appendix [X] (for 
©-independent foo(%)) in the form 



dCoo( - T) =(2z co -l)-[C oo( ,)-z ooa ;] dln(/ - (a:)) 



dx 



d.T 



(11) 



Since the solution of this differential equation is propor- 
tional to l//oo(x), we preferred to integrate Eq. (flT|) with 
the initial condition Coo(0) = / oo (0)/(l — z^) to achieve 
a stable numerical results for large x also. The resulting 
estimate for Coo(x) is shown as line in Fig. fTPTa). The 
line lies slightly above the data for the scaled capture 
numbers for T = 10 s and = 0.2, indicating that indeed 
an estimate of a limiting curve for the scaled capture 
numbers is obtained. In a cross-check, we performed the 
integral in Eq. ([7]) with the estimated (x) and recov- 
ered the line in Fig. HJa). 

For the compact island morphologies, Eq. (UJ would be 
of limited practical use, because, as discussed in Sec. [V] 
the RE ([1]) , @ fail to predict the ISD correctly already at 
small due to coalescences. Nevertheless, from a con- 
ceptual viewpoint, it is interesting to study the scaled 
ISD and their relation to the scaled capture numbers for 
the sub- islands. The corresponding data shown in Fig. II II 
indicate a behavior similar as for the fractal island mor- 
phologies, where the limiting curves are dependent on 
0. 



VII. SUMMARY 

The capture numbers entering the RE for the growth 
kinetics of thin films have been determined by KMC 
simulations in their dependence on both the coverage 
and the V = D/F ratio for the point island model 
and for two simple growth models representative for is- 
lands with compact and fractal shapes. It was shown 
that the ©-dependence of the capture numbers could not 
be accounted for by the ratio R s /£, of the mean island 
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FIG. 11. (a) Scaled island size distributions and (b) scaled 
capture numbers of sub-islands as function of the scaled island 
size in the case of compact island morphologies for fixed F = 
10 s and four coverages 0. The inset shows the corresponding 
data at fixed O — 0.2 and the four simulated F values. 



radius R s and the effective adatom capture length £ of 
the RE. This suggests that the strong deviations between 
the capture numbers determined from the simulations 
and the ones predicted by the self-consistent theory have 
their origin in the linearization step used in this theory. 
The RE with self-consistent capture numbers neverthe- 
less provide a good quantitative account of the adatom 
and island density. The deviations to the correct capture 
numbers lead, however, to a failure for a description of 
the ISD. 

Integration of the RE with the simulated capture num- 
bers determined from the KMC simulations gives an ex- 
cellent quantitative prediction of the ISDs. For the com- 
pact islands morphologies, it was found that coalescence 
events, not considered in the RE, become relevant already 
at small coverages well below ~ 0.15, where coales- 
cence events do not significantly affect the island den- 
sity. Compared to the fractal island morphologies, the 
coalescence rate for the compact morphologies is much 
higher. The ISD is affected already by a rather small 
number of coalescences, because these lead to a reshuf- 
fling of weights for different island sizes. The lower co- 
alescence rate for fractal morphologies is caused by the 
fact that fingers of two approaching fractal islands typi- 



cally first avoid each other, which subsequently leads to 
a screening effect and a slowing down of further growth 
of these fingers. 

Finally we discussed the limiting curves for the scaled 
ISDs when r — > oo. For the point islands the KMC data 
provide evidence that these limiting curves are indepen- 
dent of the coverage, which is given by the RE prediction 
([7]). This means that there exists a true scaling behav- 
ior in the T — > oo limit, where the dependence on 
is fully accounted for by the mean island size s. For 
the growth models representing compact and fractal is- 
land morphologies, the results indicate that the limit- 
ing curves are dependent on 0. This implies that one 
needs to solve the partial differential equation (|A6j) [or 
Eq. (|A10j) ] to calculate foo(x,Q) from Coa(x,&). Unfor- 
tunately no successful theory exists so far to predict the 
limiting curve Coo(x, 0) for the scaled capture numbers. 

The limiting curves are also different for different mor- 
phologies. Considering how sensitive the shape of the 
limiting curves depends on the nonlinear behavior of the 
scaled capture numbers as a function of the scaled is- 
land size, it is well possible that the shape will also vary 
with details of the growth mechanisms, even if the type 
of island morphology remains essentially the same. 



Appendix A: Rate equation prediction of the 
limiting curves for the scaled island size distribution 

For large r, s ~ A^ 1 ~ T 1 / 3 and x = s/s becomes 
a continuous variable, which allows one to derive a de- 
termining equation for the scaled ISD in dependence of 
the scaled capture numbers. The derivation was first pre- 
sented by Bartelt and Evans . 16 ' 17 Replacing the variable 
s by x and using d/d{Ft) = (1 - 0)9/90, Eqs. © can 
be written in the continuum limit as 



9? 



1 



90 



(1 - 0)5 



(Al) 



Defining f(x, 0,T) = s 2 n xS /Q, C{x, 0,T) = <j xS /<t, and 
K(x, 0, F) = Kxg/R, one has 



9 ( 

dx \ ^ xs 

dn r 



90 



e<t d(Cf) 

s 2 dx 
QR d(Kf) 

s 2 dx 
-(2*-l)/- 



(A2) 
(A3) 



where z(Q,T) = 9 In s/9 In 0. The reduced RE moreover 
predict m ~ (1 - 0)/(IW) - (1 - &)s/{@Ta) for large 
r and fixed > 0^ ~ T -1 / 2 . Inserting this relation and 
Eqs. dH El ES into Eq. flM} gives 



(2z- 



l)f + zx^-Q^l = ^Il 
dx 90 dx 



d(Kf) 



(1 - 0)5 



dx 
(A5) 
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Introducing the limits Coo(x, 9) = limr^oo C(x, 9, T), 
KagfaQ) = limr^oo 9,T) and Zoo (6) = 
limr-j-oo ^(©^r), Eq. (|A5|) yields a determining equation 
for foo(x, 9) = limr->oo f(x, 9, T). 
For limr-s-oo H/s — one obtains 



9/c 



e 



(A6) 



<9x 99 dx 

The condition limr->oo «/s = is valid for point islands, 
and it can be expected to hold also for compact island 
morphologies unless atoms deposited on top of islands 
are essentially all attaching to the island edge in the first 
layer (a situation unlikely due to second layer nucleation 
on larger islands). 

When integrating Eq. (|A6I) over x from zero to infin- 
ity, the first, second and third term on the left hand 
side yield {2z x — 1), Zoo (after a partial integration) and 
zero, respectively, because of the normalization of /oo- 
The right hand side becomes [— Coo(0, Q)f ca {Q, 9)] (note 
that for large x, C M ~ x and must decrease faster 
than x to be normalizable - simulation results show that 
foo should in fact decay much faster). Accordingly, the 
relation 



/oo(o,e) = 



i 



;(©) 



Coo (0,9) 



(A7) 



must be fulfilled. A corresponding relation can be de- 
rived in the same way already from Eq. (|A5I) . Analo- 
gously, when first multiplying Eq. (| A6|) with x and then 
integrating, one obtains 



Oo(x,9)/oo(x,9) = 1. 



(A8) 



Integrating Eq. 

Coo{x, 9) = 2 



6]) to a finite value x then yields 

l-Zoo(O) 



9 



foo(x,e) 



dz7oo(z',e) 



d 



dx' foo(x',<d) 



(A9) 



which expresses Coo{x, 9) as a functional of foo(x, 9). 

When one further assumes that the limiting curve /oo is 
independent of 9, one has dfoo/dQ = and can neglect 
the corresponding term in Eq. (|A5|) . For self-consistency, 
this requires also Coo and Zoo to become independent 
of 9. In fact, one can conversely show that if Coo and 
Zoo are independent of 9, /oo must by independent of 9 
also. Under this assumption Eq. (|A6[) then reduces to a 
separable ordinary differential equation, whose solution 
is given by Eq. ([7]), with Ctot{x) equal to Coo (a;) and z 
equal to Zoo- 

If there exists a finite limit Poo(9) = limr^oo k/s > , 
as it may be the case for fractal island morphologies (see 
the discussion in Sec. IVIj) . Eq. (|A5j) yields 



(2Zoo - l)/oo + Z 



dx 

9(Coo/oo) 

dx 



9 



99 
9 



d(K c 



(A10) 



(1-9) 



dx 



as determining equation for foo{x, 9). Strictly speaking, 
a 9-independent foo(x) should not exist then and one 
needs to solve the semi-linear partial differential equa- 
tion (|A10p . If one nevertheless makes the approximation 
dfx /dO ~ in Eq. (IA10[) and considers Coo and z x to 
be independent of (or only weakly dependent on) 9, one 
would obtain the weakly 9-dependent solution Eq. ([7]) 
with Ctot = Coo + P^QK^/il - 9). 



Appendix B: Quantitative analysis of coalescence 
events 

For a quantitative analysis of the coalescence behavior 
we determined the fraction of pair distance vectors of 
coalescing islands that before coalescence exhibit an anti- 
parallel orientation to the vector connecting the center of 
masses of the islands. Let us denote by Rjj- the vector 
pointing from the center of mass of island i to the center 
of mass of island j, and by r ia jp the vector pointing from 
atom a of island i to atom /3 of island j. The fraction of 
distance vectors with anti-parallel orientation then is 



(Bl) 



a,0 



where H(.) is the Heaviside jump function with H(x) = 1 
for x > and zero else. For a given time lag At be- 
fore coalescence, the <i>y were averaged over all coales- 
cence events, yielding the mean fraction $(At) of dis- 
tance vectors with anti-parallel orientation. To obtain 
the corresponding data, configurations generated by the 
KMC simulations were analyzed afterwards back in time, 
starting from the instant where islands first touched each 
other. 

The mean fraction $ obtained from this analysis is 
shown in Fig. [TST a) as a function of At for T = 10 7 . We 
assigned negative values to At to emphasize that 3? (At) 
was determined for lags before a coalescence event. That 
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FIG. 12. (a) Mean fraction $ of distance vectors with anti- 
parallel orientation (with respect to the center of mass dis- 
tance vector), and (b) mean minimal distance d between is- 
lands as function of the time lag At < before a coales- 
cence event at zero time. The times are given in units of F~ 
and the data were determined from the KMC simulations for 
T = 10 7 . 
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Q(At) for fractal islands is by many orders of magni- 
tude larger than for compact islands demonstrates the 
partial inter-penetration of the fractal islands before co- 
alescence. The value $(Ai) ~ 0.01 reached for the fractal 
island morphologies in the limit At — > means that on 
average about 10% of the atoms of each island in a coa- 
lescence event pass each other. That the partial inter- 



penetration is accompanied by a slowing down of the 
approach of two islands before coalescence can be seen 
in Figure I12f b). where the averaged minimal distance 
d(At) between coalescing islands is shown, that means 
dij = min aj ^(|rj a) j^|) averaged over all coalescences of 
islands i and j for time lag At. The (negative) slope of 
d(At) is significantly smaller for the fractal island mor- 
phologies, giving evidence for the screening effect fi*2I 
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